Quantum Nondemolition Charge Measurement of a Josephson Qubit 
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In a qubit system, the measurement operator does not necessarily commute with the qubit Hamil- 
tonian, so that the readout process demolishes (mixes) the qubit energy eigenstates. The readout 
time is therefore limited by such a mixing time and its fidelity will be reduced. A quantum nondemo- 
lition readout scheme is proposed in which the charge of a flux qubit is measured. The measurement 
operator is shown to commute with the qubit Hamiltonian in the reduced two-level Hilbert space, 
even though the Ifamiltonian contains non-commuting charge and flux terms. 
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Reading out the final state of qubits is an essential part 
of any quantum computation scheme. Most commonly, 
the readout is considered to be projective, meaning that 
the state of the qubit is projected onto the eigenstates 
of the measurement operator. This, however, requires 
strong qubit-detector coupling and fast readout response, 
which could be practically difficult. Furthermore, strong 
interaction with the detector may limit the coherence 
time of the system during quantum operation. For a pro- 
jective measurement, the dephasing time Tip due to the 
readout backaction should be much smaller than the free 
evolution time of the qubili(fi = 1): <C A~^, where A 
is the qubit's energy splitting. In general, however, it is 
also possible to measure the qubit's state with a weaker 
couphngSi (r^ > A"^) as long as < Tmcas < r^i^. 
The measurement time Tmcas determines when the two 
states have distinguishable output, and the mixing time 
Tmix shows when the information about the initial state 
is destroyed due to the interaction with the detector. In 
this regime (called "Hamiltonian dominated" in Ref. 0), 
the qubit is measured in the eigenstates of its own Hamil- 
tonian and the detector's output is the expectation value 
of the measurement operator in these eigenstates. The 
backaction noise causes extra relaxation of the qubit to 
the ground state, limiting Tmix. The reported single shot 
readout of fiux qubits using DC- SQUID in Ref. ^ seems 
to agree more with this picture than with the projective 
measurement. 

In general, it is desirable to have a long Tmix to allow 
a long measurement time with high readout fidelity. In a 
Quantum Nondemolition (QND) measurement, while the 
readout process collapses the state of the qubit, it does 
not demolish it.^ In other words, the readout measures 
the eigenstates of the qubit Hamiltonian without chang- 
ing them (e.g. causing fast relaxation). Mathematically, 
it means that the measurement operator commutes with 
the qubit Hamiltonian. Hence, Tmix will be limited only 
by the intrinsic relaxation time of the qubit, which could 
be long. While increasing the readout fidelity, QND mea- 
surement allows reducing the qubit-readout coupling and 
therefore decreasing unwanted decoherence due to inter- 
action with the detector during the operation time. 

Despite the importance of QND measurement, there 
are only few QND readout schemes for superconduct- 
ing qubitsi^ In the Saclay's charge-flux qubit design 
(quantronium) jfi the measurement operator (flux) com- 
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FIG. 1: Three- Josephson-junction charge-flux qubit coupled 
to an SET for QND readout. 



mutes with the qubit Hamiltonian only at the flux de- 
generacy point. At the time of readout, however, the 
qubit is moved away from this point, hence resulting in a 
readout that is not a QND type. This may be responsi- 
ble for the observed small readout fidelity. Another type 
of charge-flux qubit, which is flux based, was proposed 
in Ref. T". Like in the quantronium case, the qubit can 
be optimally protected against dephasing due to low fre- 
quency noises at the charge-flux double degeneracy (so 
called magic) point. In this paper, we suggest a modifi- 
cation to the qubit, to make QND readout possible. 

The qubit (Fig. consists of a superconducting loop 
containing three Josephson junctionSj-S' threaded by an 
external flux close to half a flux quantum ^ ^o/2 — 
h/Ae). The Josephson energy and junction capacitance 
of two of the junctions are the same {Eji = Ej2 = Ej, 
Ci — C2 — C) while those of the third junction are 
slightly smaller [aEj and aC, with 0.5 < a < 1). In ad- 
dition, two voltage sources VgA and VgB are capacitively 
connected to two of the islands, which are also coupled 
capacitively to an SET or RF-SET— as a sensitive charge 
detector for readout. At an appropriate biasing point, 
state dependent voltages appear on these islands^i which 
induce a charge on the SET's island, affecting the current 
through it. 

We define (f) = -|-(/)2)/2 and 6* = (0i — (/'2)/2, where 
01,2 are the phase differences across the two larger junc- 
tions. The Hamiltonian of the system is H = Hq + 
HsET + ^^int , where i?sET is the SET Hamiltonian, and 
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is the qubit Hamiltonian. Here 

Ui(j),9) = Ej[-2cos(j)cose + acos{2TTf + 29)] 

is the potential energy, — ~id/d^ and Pg = —id/dg 
are the momenta conjugate to (j) and 9, — C(l + 
7)/2e2, Mg = C(l + 7 + 2a)/2e2, 7 = C^/C, and 
n± = ± "-gs, with UgA^gB = VgA,gBCg/2e being the 
normahzed gate charges. Throughout this paper we take 
/ = $a;/$o - 1/2 = 0, leading to a U{(j),9) with degen- 
erate minima at (/> = 0, = ± arccos(l/2a). 
The interaction Hamiltonian is written as 



(1) 



where k± = {CcA i Ccb)/Cy., n represents the normal- 
ized charge (operator) of the SET island, and Cs is its 
effective capacitance. 

The commutation relation 

[Hiat.HQ] = -i(sin(/)C0s6l)-^n 

— ifcosoi<sin6' — 2Q!sin2^?) h 

^ 'Mg 

is clearly nonzero, meaning that in general the readout 
can change the eigenstates of the qubit Hamiltonian. For 
a qubit well in the flux regime, however, the qubit wave 
function is confined within a unit cell in the phase space. 
The fluctuations of </> is therefore small and, for the low- 
est energy eigenstates, w 0. Thus, [i7int,i?Q] oc k_ 
vanishes if CgA — CgB, suggesting a scheme for QND 
measurement. However, inter-unit-cell tunnelling, which 
is important for operation in the charge-flux regime^ will 
delocalize 0, conflicting with the above picture. There- 
fore, it is not possible to have -ffint exactly commute with 
Hq. As we shall see however, it is possible to make them 
commute in the Hilbert space reduced to the first two en- 
ergy levels, hence satisfying a much weaker requirement: 
(l|//int|0) = 0, yet sufficient for QND measurement. The 
backaction noise in that case will not affect the rate of 
relaxation between the two lowest energy states (|0) and 
|1)) of the qubit during readout, i.e. will not reduce r^iy^. 

To study in more detail, we calculate the relaxation 
rate due to the SET's backaction on the qubit: 
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[5(A) + 5(-A)], (2) 



where PVg = {i\P^^g\j), and 



S{uj) = / e''^\n{t)n{Q))dt (3) 



is the noise spectral density. Using the orthodox theory 
for a normal SETjii the spectral density S'(w) at small 
bias current is given bjii^ 
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where / = erir2/(ri-|-r2) is the average current through 
the SET and Fi (F2) is the rate of tunnelling of an elec- 
tron into (out of) the island via the first (second) SET 
j miction: 
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ASi,2 = ±2Ec^[ngs- 1/2) + [l/1)eVDS- (6) 

Here R is the normal resistance of the junctions, -Ec^ = 
e^/2CE is the charging energy, ng^ = CgsVgs/e is the 
normalized gate charge, and Vds is the voltage between 
drain and source of the SET. 

Let us first consider the case where the readout is cou- 
pled to only one of the islands (as in Ref. 0): CcB — 0, 
therefore = k_ = Cca/Cy.- An SET, operating in 
the quasiparticle branch, typically has 7^1 nA and 
Vds - 400 //V, which lead to Fi,2 lO^'' s^^. For 
a typical qubit with Eq/Ej = 0.1 and a — 0.75, as 
suggested in Ref. H we find: « 5 GHz, A « 7 

GHz, « 2.5M-1 w AEc w 20 GHz, P^^o « 0.04, 

and Pg" « 0.5. Substituting into Q and ^ and using 
Cca/Cy, ~ 0.2, we find a mixing time of Tmix ^ 20 ns, 
which is probably too short to read out the qubit. It is 
therefore important to use symmetric coupling (see be- 
low) or operate the SET at the double Josephson quasi- 
particle point, "'^^ where the SET's backaction noise (but 
also its sensitivity) is significantly smaller.^" Notice that 
Pg° ^ P^^ (see also Fig. |21l, hence the second term in 
has much larger contribution to the backaction noise 
than the first one. 

In the symmetric coupling case {CcA — Ccb), — 
and therefore the second term in Q is cancelled. The 
relaxation rate will then only depend on P^'', which is 
typically much smaller than Pg^. As we shall see, there 
even exist special points at which P^° = 0, hence Q 
vanishes, making QND readout possible. Let us write 
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d9-^i{^,9){~id^)-i>o{^,9), (7) 



where ^0,1 are the flrst two energy eigenfunctions. When 
n-gA = rigB — 0, the wave functions have the following 
symmetries 



*o,i(-0,^) = ^o,ii(t>,0) 
^oA^,-0) = ±*o,i(0,e) 



(8) 
(9) 



Either of these is enough to make ^ vanish. The output 
however will be zero at this "magic" point. It is therefore 
necessary to move the qubit away form this point at the 
time of readout. As we shall see, unlike in the quatronium 
case, it is possible to flnd points at which O vanishes 
while the measurement output is finite. 

To see explicitly how the symmetries © and © 
change at finite UgA and Ugs, we first consider the sim- 
pler problem of a one dimensional particle in a symmetric 
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FIG. 2: \Pl° \ and \Pi°\, plotted vs UgA when Ugs = 0.25 (a) 
and UgB = 0.5 (b), for a qubit with Ec/Ej = 0.2 at / — 0. 
The top curves are for \Pe'^\ and the bottom ones show l-P^"]. 
The values of a are 0.6 (solid line), 0.7 (deshed line), 0.8 
(dotted line), and 0.9 (dash-dot line). 

periodic potential V{x) 

H=^{P + eA/cf + V{x), (10) 
Zm 

V{x + d) = V{x) = V{~x), (11) 

where A is a uniform vector potential. The eigenfunc- 
tions of this system are described by the Bloch wave 
functions. The ground state of the system has zero crys- 
tal momentum and is therefore the solution of (|10|l with 
periodic boundary condition \l/o(a; + d) — 'i'o{x). When 
A = 0, the ground state is symmetric: ^'o(— a;) = ^"0(2;). 
At finite A, one can perform a gauge transformation to re- 
move A from H1U|I . The boundary condition then becomes 
^Q{x + d) = e'x^'o(x), where x = j{e/hc)Adx = eAd/hc. 
When the wave function is confined within the wells, a 
standard tight binding type of calculation can be used to 
show that a small antisymmetric term, proportional to 
sinx, will be added to the otherwise symmetric ground 
state. In general, the wave function of the system in the 
j-th state, has the form 

^,ix)^^°ix) + ^l{x)smx, (12) 

where ^l{x) is a small term with opposite symmetry 
as that of the main term \I'"(x). '^}{x) is dominantly 
determined by the geometry of the well and therefore 
is independent of x- For example, when the two low- 
est levels of the system are well separated from the rest, 
^l{x) cx «'?(a;). 

Let us now return to our qubit system. The two di- 
mensional potential energy of the system is 27r periodic 
in both and (j)2 directions. In terms of (j) and 9, the 
boundary condition reads 

*(</> + vr, 6* ± tt) = e"("+±"-)*((/., 9). (13) 

A straightforward generalization of the above argument 
yields 

^,{^,9)=^°i(l),9)+^licl),9) sin7rn+ 

+'i'f{(f>, 9) sin7rn_ -I- (0, 9) sin7rn_|_ sin7rn_. 



For the ground (first excited) state [i = (1)], and Vtf 
are symmetric (symmetric) functions of (j), while ^\ and 
'^^ are antisymmetric (antisymmetric). Likewise, and 
'^l are symmetric (antisymmetric) in 9, while and ^'f 
are antisymmetric (symmetric). The difference between 
(j) and 9 stems from the fact that the potential has double- 
well structure only in 9 direction. Neglecting the small 
contribution from the decay region, equation Q gives 

P^^ — Zsin7m+ sin7m_ 

= -Z(cos 27rng^ — cos 27rngB), (14) 

where Z ^ All + A\l + A\\ + Ajg, and 

'^'t'j d9^n<l^.0){-id^)^^^{cp,9). 

Similarly, for Pg we find 

Pgio = + B\l sin^ 7m+ + BH sir? 7rn_ 

-I-B33 sin^ 7™+ sin^ 7rn_ , (15) 

where 

'^'t'j d9^fi^,9)i-zde)^^i<p,9). 

While Pg" is always finite, P^" vanishes when either n+ 
or n_ is an integer number. At these special points, 
(l|7?int|0) = if /t_ = in Ql. The latter can be easily 
achieved by symmetrically coupling the two islands of the 
qubit to the readout SET {CgA — Cgs)- 

The above argument is valid beyond the approxima- 
tions used to derive (fT^ and (fT^ (see Fig. [5)). Indeed, 
symmetries ^ and © always hold when respectively 
and n_ are integer numbers, causing (|7|) to vanish. A 
simple way to access these symmetry points is to adjust 
the gate voltages so that ngA = ±ngB, which correspond 
to rizp = 0, respectively. 

Figure |21 shows the result of numerical calculation of 
I ^0*9 1 fo'' different values of a. When a is small, the wave 
function of the qubit is localized within a unit cell. As 
a result, the jP^^I curve agrees very well with (|14|l . A 
larger a delocalizes the states, causing deviation from a 
simple sinusoidal behavior. When UgB = 0.25 (Fig. |2^), 
jP^"! vanishes at UgA — 0-25 (n_ — 0) and UgA = 0.75 
(n+ — 1), in agreement with the above argument. For 
ngB = 0.5 (Fig. Els) both points coincide at UgA — 0-5 
(n_ — and — 1). 

While the points UgA — ngB and UgA = —ngB both 
result in P^*^ = 0, only the former is suitable for readout. 
To see this, notice that the detector measures the expec- 
tation value of the sum of the island voltages in the |0, 1) 
states: (dH/dn+) = Pf/M^, i = 0, 1. However 

P«° ^ sin.n+[A^; + ATo + {A^i + A^§) sin^ nn^], 
P" = simrn+iA'al + A^ + (A^ + A^) sin^ 
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FIG. 3: P" - P°° vs n+ for the qubit system of Fig.Elwhen 
n_ — 0. Different lines correspond to different a the way 
described in the caption of Fig. |21 

both vanish when = {ugA = —UgB), resulting in 
indistinguishable (zero) outputs. Physically, the voltages 
of the two island will be exactly equal but with opposite 
signs, cancelling each other's effects on the output. Thus, 
the qubit should be read out at n_ — (or any integer 
number). The difference in the output between the two 
states is then proportional to 

- = sinnn+iAll + A\l ~ ^[J? - ^?[J], (16) 

which vanishes at integer n-|_, as expected. 

Figure O displays the result of numerical calculation 
of — for different values of a, for the same set 
of parameters as in Fig. [2| The plot is vs n+ while 
n_= {rig A = TigB)- The sinusoidal- type behavior, is 
in agreement with (|15|l . Maximum output is achieved 
when rij^K, 0.5 {rig a = rigB ~ 0.25). This is indeed the 
optimal point for reading out the qubit. 

None of the above arguments depended on the qubit 's 
regime of operation, although we mainly focussed on the 
flux regime. The only restriction is that the qubit is- 
land voltages should not be too small to be detectable 
by the SET; the qubit should not be too far in the flux 



regime. It is indeed interesting to see how the scheme 
works when the qubit is in charge regime {Eq Ej). 
Let us write the qubit wave function in charge basis 
InA^ns), where ua.b is the number of Cooper pairs 
on the two islands. When UgA = rigB {ri- = 0), the 
states |n, n -I- 1) and \n + will be degenerate. The 
presence of the Josephson term of the third junction 
will remove this degeneracy, making the rotated states 
|n,±) = (|n, n+ 1) ± |n + 1, n))/-\/2 preferred basis for the 
qubit. In the Hamiltonian matrix, |n, has nonzero off- 
diagonal matrix elements with other charge states, hence 
will mix after complete diagonalization. The matrix el- 
ements of \n, —) with other states, however, are all zero 
as long as Eji = Ej2- It fact, \n, — ) is an eigenstate of 
the qubit Hamiltonian, as well as an eigenvector of Hint 
if K- = 0. Thus, Hint will not mix |n, — ) with any other 
states. It is easy to show that for all values of n+, one of 
the two lowest energy levels of the qubit will be a \n,—) 
state (with an appropriate value of n depending on n^), 
hence (l|i?int|0) = 0. 

To summarize, we have shown that QND measurement 
of a three Josephson junction qubit is possible if both is- 
land charges are measured symmetrically at a symmetric 
biasing point. The scheme is applicable for a wide range 
of parameters from charge to flux regime as long as charge 
measurement is possible. QND readout allows long read- 
out time with high fidelity and small qubit-detector cou- 
pling to prevent decoherence due to the interaction with 
the measurement device during quantum operation. It 
also makes it possible to operate the SET in the quasi- 
particle branch benefitting from larger charge sensitivity. 
Experimental study of such a system can shine light on 
the intrinsic mechanisms of relaxation and decoherence 
sources in charge or charge-flux qubits. 
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